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Abstract 

We present an analysis leading to a conjecture on the exact location of the multicritical 
point in the phase diagram of spin glasses in finite dimensions. The conjecture, in satisfac- 
tory agreement with a number of numerical results, was previously derived using an ansatz 
emerging from duality and the replica method. In the present paper we carefully examine 
the ansatz and reduce it to a hypothesis on analyticity of a function appearing in the duality 
relation. Thus the problem is now clearer than before from a mathematical point of view: 
The ansatz, somewhat arbitrarily introduced previously, has now been shown to be closely 
related to the analyticity of a well-defined function. 

1 Introduction 

To establish reliable analytical theories of spin glasses has been one of the most challenging 
problems in statistical physics for years. The problem was solved for the mean- field model 0121 
0111]. Much less is known analytically for finite-dimensional spin glasses, for which approximate 
methods including numerical simulations and phenomenological theories |2| have been the main 
tools of investigation in addition to a limited set of rigorous results and exact solutions [HI [7] . 

Our main interest in the present contribution does not lie directly in the issue of the properties 
of the spin glass phase. We instead will concentrate ourselves on the precise (and possibly exact) 
determination of the structure of phase diagram of finite-dimensional spin glasses. This problem 
is of practical importance for numerical studies since exact values of transition points greatly 
facilitate reliable estimates of critical exponents in finite-size scaling. 

More precisely, recent developments [HllHlEl^DEl to derive a conjecture on the exact loca- 
tion of the multicritical point in the phase diagram will be analyzed from a different view point. 
The conjecture was derived using the replica method, gauge symmetry and a duality relation. 
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In addition it was necessary to introduce an ansatz to identify the location of a singularity of 
the free energy. The resulting conjecture for the transition point (multicritical point) is nev- 
ertheless in satisfactory agreement with a number of numerical results, which renders a strong 
support to the validity of our prescription. In the present paper we present a more systematic 
analysis leading to the ansatz, thus reducing the problem to that of the analyticity proof of a 
function related to duality. The ansatz was introduced somewhat arbitrarily previously. The 
discussions in the present paper make it clear that the ansatz is closely related to the analyticity 
of a well-defined function, which paves a path toward the formal proof of the conjecture. 

2 Duality relation for the replicated system 

For simplicity, let us consider the ± J Ising model on the square lattice. It is possible to apply 
the same line of argument as presented below to other systems (non- Ising models and/or other 
lattices) as will be mentioned in the final section. The Hamiltonian is 

H = -jY^rijSiSj, (1) 

where Tjj = ±1 is a quenched random variable with asymmetric distribution and J > 0. Periodic 
boundary conditions are imposed. We accept the replica method in this paper and do not strive 
to rigorously justify the validity of taking the limit n ^ in the end, where n is the number of 
replicas. 

The n-replicated partition function after configurational average is a function of edge Boltz- 
mann factors: 

[Z''] = Zn{xo{K,Kp),xi{K,Kp),--- ,Xn{K,Kp)). (2) 

Here the square brackets denote the configurational average, K stands for pj = J/ksT, and Kp 
is a function of p (the probability that Tij is 1) defined through e~^^p = (1 —p)/p- The kth. edge 
Boltzmann factor Xk{K, Kp) (A; = 0, 1, • • • , n) represents the configuration-averaged Boltzmann 
factor for interacting spins with k antiparallel spin pairs among n nearest-neighbor pairs for a 
bond (edge) as illustrated in Fig. 

Xk{K, Kp) = pe("-2*^)^ + (1 - p)e-("-2*^)^. (3) 

The expression on the right-hand side of Eq. ((2j emphasizes the fact that the system proper- 
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ties are uniquely determined by the values of edge Boltzmann factors because the interactions 
(or equivalently, the edge Boltzmann factors) do not depend on the bond index (ij) after the 
configurational average. 

The formulation of duality transformation developed by Wu and Wang J3] can be applied 
to the function Z„ of Eq. ((2)) to derive a dual partition function on the dual square lattice with 
dual edge Boltzmann factors. The result is, up to a trivial constant, [S] 

)• (4) 

The dual Boltzmann factors are discrete multiple Fourier transforms of the original Boltzmann 
factors (which are simple combinations of plus and minus of the original Boltzmann factors in 
the case of Ising spins) : 

n 

= -D^ Xk, (5) 

k=0 

where the Z?^ are the coefficients of the expansion of (1 — t)"^{l + t)""™ : 

n 

{i-tr{i+tr-"^ = Y^Dit^ (6) 

A:=0 

or explicitly, 

k 

\ ^ . _ I m \ I n — m 



For example, in the case of n = 2, the dual Boltzmann factors are given as 

2x1 = (^0 + Xi) + (Xi + X2) = Xq + 2xi + X2 

2x1 = (^0 - Xi) + (Xi - X2) = Xq- X2 (8) 
2X2 = (Xo - Xi) - (Xi - X2) = Xo - 2xi + X2. 



t t 
t t 

t- i 

t t 



Figure 1: The edge Boltzmann factor Xk{K,Kp) represents the weight for k antiparallel spin 
pairs among n nearest-neighboring pairs. The case of /c = 1, n = 4 is shown here. 
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It will be useful to measure the energy from the all-parallel spin configuration {k = 0) 
and, correspondingly, factor out the principal Boltzmann factors, xq (left hand side) and Xq 
(right hand side), from both sides of Eq. @. Since the partition function is a homogeneous 
multinomial of edge Boltzmann factors {xq,xi, ■ ■ ■ ,Xn) of order Nb (the number of bonds), the 
duality relation Q can then be rewritten as 

{xo)^''Zn{ui,U2, ■■■ ,Un) = {Xq)'^'' Z„{ul, U2, • • • , <), (9) 

using the normalized edge Boltzmann factors and (A; = 1, 2, • • • , n) defined by = x^/xq 
and ul = xl/xl. 

The discussions so far have already been given in references [HI IIUI ITT] . In those papers we 
went on to try to identify the multicritical point by the fixed-point condition of the principal 
Boltzmann factor, XQ(K,Kp) = Xq{K, Kp), combined with the Nishimori line (NL) condition 
K = Kp iQ. The reason for the latter choice is that the multicritical point is expected to lie 
on the NL |14j . In this way an ansatz was made at this stage that the multicritical point is 
given by the relation xq{K,K) = Xq{K,K). The resulting expression for the location of the 
multicritical point was confirmed to be exact in the cases of n = 1, 2 and n ^ oo. Extrapolation 
to the quenched limit n ^ gave results in agreement with a number of independent numerical 
estimates as listed in Table 1. However, it was difficult to understand the mathematical origin 
of the somewhat arbitrary-looking ansatz. 

We develop an argument in the next section to justify the above-mentioned relation xq{K, K) = 
Xq{K,K) as the fixed point condition of duality relation for the replicated it J Ising model on 
the NL. 

3 Self-duality 

The duality relation © applies to an arbitrary set of parameter values {K, Kp) or an arbitrary 
point in the phase diagram (Fig. [21). However, we restrict ourselves to the NL (K = Kp) to 
investigate the location of the multicritical point. Then Eq. Q is expressed as 

{xo{K))^^Zn{ui{K),U2{K),--- ,Un{K)) = [xliK))'"'^ Zn{u\{K),ul{K), - ■ ■ ,ul{K)) (10) 

since xq, ui, • • • , and Xq, uj, • • • , u* are now functions of a single variable K. 
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It is difficult to apply directly the usual duality argument (identification of a fixed point with 
the transition point assuming uniqueness of the latter) to Eq. (|inj) . The reason is that Zn is a 
multi-variable function: The two trajectories representing C{K) *== {ui{K),U2{K), ■ ■ ■ ,Un{K)) 
and C*{K) {u{{K),U2{K), ■ ■ ■ ,u!^{K)) do not in general coincide as depicted in Fig. |31 In 
other words, there is no fixed point in the conventional sense for the present system. 

In spite of this difficulty of the apparent absence of a fixed point, we now show that it is 
still possible to devise a map which renders the present replicated system self-dual with a fixed 
point. To derive such a result, it is useful to note a few properties of the quantities appearing 
inEq. (fTn|) . 

Lemma 1 The normalized Boltzmann factor Uk{K) is a monotone decreasing function of K 
from Uk{0) = 1 to limA'^oo 'Ufc(i^) = 0. The dual u1{K) is a monotone increasing function from 
ul{0) = to limii-^oo ul{K) = 1. Here G {1, 2, • • • n} and n G N. 

Proof. From the definition (0) of and the NL condition K = Kp (i.e. e~^^ = (1 — p)/p), it 
is straightforward to verify that 

^(n+l-2k)K _^ ^-{n+l-2k)K 
^fc(^) = e(n+l)X + ^-{n+l)K ' (11) 

from which the result for Uk{K) follows. The dual are shown to satisfy [HI 



(tanhifj'^+i (/c odd) 
(tanhi^)'^ (A: even) 



(12) 



which leads to the statement on u^. 



Lemma 2 The normalized partition function Zn{ui, ■ ■ ■ ,Un) is a monotone increasing contin- 
uous function of all of its arguments ui, - ■ ■ ,Un with the limiting values in the hypercube [0, 1]"' 
Zn{0, ■ ■ ■ , 0) = 2 and Zn{l, ■ ■ ■ , 1) = 2"^, where N is the number of sites. 

Proof. Since the normalized partition function is a multinomial of normalized edge Boltzmann 
factors ui, - ■ ■ ,Un with positive coefficients, the first half of the above statement is trivial. 
To check the second half, we note that ui = • • • = n„ = corresponds to the case where no 
antiparallel spin pairs are allowed in any replica at any bond. The only allowed spin configuration 
is the all-parallel (i.e. perfectly ferromagnetic) one, for which we have set the energy (or the 
edge Boltzmann factor uq = 1) since the energy is measured from such a state by dividing by 



xq. Taking into account the global inversion degeneracy, we conclude that -E„(0, • ■ ■ ,0) is equal 
to 2uo = 2. 

Similarly, when ui = ■ ■ ■ = Un = 1 (corresponding to the high-temperature limit), all spin 
configurations show up with equal probability. Therefore the normalized partition function just 
counts the number of possible spin configurations, yielding 2"^^ . | 

Now we are ready to step toward the main theorems. 

Theorem 1 There exists a monotone decreasing function d{K) with \miK—>Q d{K) oo and 
limK^ood{K) = by which the value of Z at {ul{K), ■ ■ ■ ,U^{K)) becomes equal to the value at 

{ui{d{K)),--- ,Un{d{K))), 

ZM{K),--- ,uUK)) = Zn{ui{d{K)),--- ,Un{d{K))). (13) 

Proof. According to Lemma the curve C{K) = {ui{K), ■ ■ ■ ,Un{K)) starts from the point 
(1, • • • ,1) and ends at (0, • • • ,0) as K increases from to oo. Similarly, the curve C*{K) = 
{u\{K),--- , ^^^(/C)) starts from (0, •• • , 0) and ends at (1, • • • ,1). Thus ^„(iii(K), • • • ,n„(i^))'= 
K{K) is continuous and monotone decreasing from Z//(0) = 2"^ to \\m.K^coU{K) = 2 by 
Lemma 121 Similarly, Zn{u\{K),--- ,u'^{K)) *== V{K) is continuous and monotone increasing 
from V(0) = 2 to limx->oo V(i^) = 2"^. Consequently there exists a monotone decreasing func- 
tion d{K) (satisfying limi^-^o d{K) oo and limK^ood{K) = 0) that relates the values of lA 
and V such that U{d{K)) = V{K). I 

Corollary 1 The normalized partition function satisfies the duality relation 

{xo{K)f^Zn{ui{K),--- ,Un{K)) = {xl{K)fBZn{ui{d{K)),--- ,Un{d{K))), (14) 

where d{K) is the monotone decreasing function shown to exist in Theorem 1. 

Proof. Immediate from Eq. (|1U() and Theorem 1. 

Theorem 2 The normalized partition function is self-dual with a fixed point Kq given by Kq = 
d{Ko), which is equivalent to Xo{Kq) = Xq{Kq). 

Proof. The first half of the statement is immediate from Corollary 1. The second half comes 
from the observation that the values of Zn on both sides of Eq. ()14() become identical at the 
fixed point Kq = d^Ko), thus making the prefactors, xo{K) and Xq(K), equal to each other. 
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Corollary 2 Assume that d[K) has no singularity for K G [0, oo). // the free energy per site 
of the replicated system has a unique singularity at some K = in the thermodynamic limit, 
then Kc is equal to Kq of Theorem 2. 

Note. As long as d[K) is analytic, the above statement is the same one as the usual duality argu- 
ment for the ferromagnetic Ising model on the square lattice, in which d{K) is —(1/2) log tanh K. 
However, if d{K) happens to be singular at Ki for example, this singularity would be reflected 
in the singularity of the free energy (away from the fixed point Kq) through the relation (|14() : 
The free energy per spin derived from the right-hand side will be singular at Ki reflecting the 
singularity of d{K) there. This causes a singularity of the free energy derived from the left-hand 
side at Ki. 

We have been unable to prove analyticity of d{K). This is one of the reasons that the 
following two statements are conjectures. The other reasons include the validity of the replica 
method and the absence of a formal proof for the existence of the multicritical point on the NL. 

Conjecture 1 [HI E] The exact location of the multicritical point for the n-replicated it J Ising 
model on the square lattice is given by the relation xq = Xq with K = Kp, i.e. 

Note. The explicit expressions of xq and Xq are given in Eqs. (jSJ and ©. See also references 

puniini. 

Conjecture 2 [UllHl The exact location of the multicritical point for the it J Ising model on the 
square lattice with quenched randomness is given by the formula 

-Pclogpc - (1 - Pc) log{l - Pc) = ^^^. (16) 
Proof. The limit n ^ of Eq. (|15|) and the NL condition K = Kp yields the above formula. 

4 Concluding remarks 

The main physical conclusions. Conjecture 1 and Conjecture 2, are not new. The significance 
of the present paper is that we have derived the ansatz to identify the multicritical point, 
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xq{K) = Xq{K), from the assumption of analyticity of the function d{K), thus hopefully coming 
a little closer to the formal proof of the conjecture. 

In this paper we have limited ourselves to the itJ Ising model on the square lattice for 
simplicity. It is straightforward to apply the same type of argument to other lattices and other 
models. For example, models on the square lattice (such as the Gaussian Ising spin glass and the 
random chiral Potts model) can be treated very similarly: The differences lie only in the explicit 
expressions oi as given in section 2.9 of reference ,9, and section 4.1 of reference Also, the 
duality structure of the four-dimensional random plaquette gauge model is exactly the same as 
the lb J Ising model on the square lattice ^0], and therefore the present analysis applies without 
change. In the case of mutually dual pairs of lattices, such as the triangular and hexagonal 
lattices or the three-dimensional it J Ising model and the three-dimensional random-plaquette 
gauge model, a simple generalization suffices that refers to the composite duality relation of 
the two systems, xo{Ki)xo{K2) = Xq{K2)xq{Ki), where Ki is the critical point for one of the 
systems and K2 is for its dual, as detailed around Eq. (17) of reference |llj . 

Let us next give a few remarks on the comparison with numerics in Table 1. As is seen there, 
our conjecture agrees with a number of numerical results but lies slightly outside error bars for 
some instances for the square lattice, 0.8907(2) 17^ and 0.8906(2) W vs. 0.889972 of our 
conjecture. A similar situation is observed for three pairs of mutually dual hierarchical lattices 
analyzed in j25j : The sum of the values of binary entropy H{p) = —p log2 p— (1 —p) log2(l — p) for 
the pair of mutually dual multicritical points is exactly equal to 1 according to our conjecture, 
whereas the numerical results are not precisely unity, 1.0172, 0.9829 and 0.9911. We have 
no definite ideas at the present moment where these subtle differences come from. Further 
investigations are necessary. 

The final remark is on the transition points away from the NL (or the shape of the phase 
boundary away from the multicritical point). We expect (but cannot prove) that the rela- 
tion xo{K,Kp) = Xq{K,Kp) gives the true critical point only on the NL. An important rea- 
son is that the limiting behavior of the transition point as p ^ 1, derived from the ansatz 
xo{K,Kp) = Xq{K,Kp), shows a deviation from a perturbational result, see 0. The particu- 
larly high symmetry of the system on the NL ^2Xi could be a reason for the success only on the 
NL. 
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Figure 2: A typical phase diagram of the it J Ising model. Shown dash-dotted is the NL. Note 
that the spin glass phase is believed not to exist in two dimensions. The multicritical point is 
denoted as MCP. 




Figure 3: Trajectories of the points C[K) = {ui{K),U2[K)^- ■ ■ ,Un{K)) and C*{K) = 
{u\{K),U2[K)^- ■ ■ ,u'^{K)), projected onto the {ui,U2) plane. The arrows indicate the sense 
of motion of C{K) and C*{K) as K changes from to oo. 
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Model 


Numerical estimates 


Reference 


Our conjecture 




Reference 


SQ Ising 


0.8900(5) 




m 


0.889972 




[Hlil 




0.8894(9) 




m 










0.8907(2) 














0.8906(2) 




CHI 










0.8905(5) 




m 








SQ Gaussian 


1.00(2) 




m 


1.021770 




[HI El 


SQ 3-Potts 


0.079-0.080 




m 


0.079731 




[HI El 


Ad gauge (RPGM) 


0.890(2) 




m 


0.889972 




uni 


TR 


0.8355(5) 




m 


0.835806 




m 


HEX 


0.9325(5) 




m 


0.932704 




m 


3d Ising (RBIM) 


0.7673(3) (= Pa) 




121 


— 






3d gauge (RPGM) 


0.967(4)(= p,2) 




m 


— 






RBIM+RPGM 




= 0.99(2) 




H{pci) + H{p,2) 


= 1 


m 


Hierarchical 1 (HI) 


0.8265(= Pea) 




m 


— 






Dual of HI (dHl) 


0. 93380 (= Pca) 




m 


— 






Hl+dHl 


H{pc3) + H{pc4) 


= 1.0172 




Hipcs) + H{pc4) 


= 1 


m 


Hierarchical 2 (H2) 


0.8149(= pc5) 




m 


— 






Dual of H2 (dH2) 


0.94872(= p,e) 




m 


— 






H2+dH2 


H{pc5) + Hip^e) 


= 0.9829 




H{pc5) + H{pc6) 


= 1 


m 


Hierarchical 3 (H3) 


0.7527(= pc7) 




m 








Dual of H3 (dH3) 


0.97204(= p,s) 




m 








H3+dH3 


H{p,7) + H{p,s) 


= 0.9911 




H{pc7) + H{p,s) 


= 1 


m 


Hierarchical 4 


0.8902(4) 






0.889972 




jsnni 



Table 1: Location of the multicritical point by recent numerical studies and our conjecture. SQ 
stands for the square lattice, and TR/HEX for the triangular /hexagonal lattices, respectively. 
RBIM stands for the random-bond Ising model and RPGM is for the random-plaquette gauge 
model. See reference ^5J for the explicit definition of three types of hierarchical lattices (HI, 
H2, H3) and their duals. Hierarchical 4 is a self-dual lattice, for which the analysis in ^1 
applies directly (though not stated explicitly in these references). The values are for pc of the 
=b J model except for the Gaussian randomness on the square lattice for which the values 1.00(2) 
and 1.021770 are for Jqc/J- Spin variables ar43sing excepting the three-state Potts model as 
indicated. The symbol H{p) is for the binary entropy H{p) = — plog2P — {1 — p) log2(l — p)- 



